Jafarkhani-Tarokh 2000-0269 

Multiple Transmit Antenna Differential Detection from 
Generalized Orthogonal Designs 

Related Applications 

5 This is a cootinuaticm-ia-part application of an application that bears the serial 

number 09/296,514, which was filed on April 22, 1999, and which is incorporated by 
reference herein This application also claims priority from Provisional Application 
60/201,342, filed May 2, 2000. 

10 Background 

This invention relates to space- time coding and, more particularly, it relates to 
space- time coding and transmission in an environment that has no knowledge of the 
transmission channel parameters. 

Hocwald and Marzetta, in "Unitary space- time modulation for multiple- antenna 

15 commvuiications in Rayleigh flat fading," IEEE Trans. Inform. Theory, March 2000, 
proposed unitary space-time codes with interesting codes, but they have exponential 
encoding and decoding complexity. Subsequently, they came up with a second 
construction that has polynomial encoding but exponential decoding complexity, which 
makes their use formidable in practical situations. In the aforementioned 09/296,5 14 

20 patent application, we proposed a coding scheme based on orthogonal designs, it being 
the first scheme that provides simple encoding^decoding algorithms. Thereafter, Hughes 
introduced a construction based on group codes that also has simple encoding/decoding 
algorithms, and Hocwald, et al. presented a construction with polynomial encoding but 
with a suboptimal decoding algorithm. 

25 This disclosure is provided because the parent application discloses an 

arrangement where two transmitting antennas are employed and at least to some artisans 
questioned whether the arrangement can be expanded to more than two antennas. 

Summary 

30 A differential detection approach allows communication without knowing the 

channel parameters, where the transmitter can include a plurality of transmitting antennas 
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ill excess of two. Each block of bk bit is mapped to a vector and the vector is processed 
to develop a set of symbols for the block, by employing a set of mutually orthogonal 
vectors and symbols developed for the immediately previous block. The set of symbols 
is then applied to a space- time mapping block, which develops signals that are applied to 
5 the plurality of transmitting antennas. 

At the receiver, for each receive antenna the received signals of a block are 
formed into vector. That vector, and a corresponding vector of the immediately previous 
block are processed in accordance with the mutually orthogonal vectors used by the 
receiver to create a combined vector. The combined vector is summed with combined 
10 vectors developed for the received signals of the other receive antennas, and the summed 
result is applied to a minimum distance decoder module. The decoder module selects a 
"most likely" vector, and that vector is applied to a mapping that is die inverse of the 
13 mapping used at fee input of the transmit 

15 Brief description of the Drawing 

}j FIG. 1 presents a block diagram a transmitter in accord with the principles of this 

i : i. invention; and 

I ^ FIG. 2 presents a block diagram of a receiver in accord with the principles of this 

*P invention 

I j 20 

M Detailed Description 

The following considers a multiple antenna, wireless, communication system like 
the one depicted in FIG. 1, under the assumption that fading is quasi-static and flat, where 
a base station 10 is equipped with n transmitting antennas and a remote receiver is 

25 equipped with m receiving antennas. At each time slot t , signals c t , i = 1 ,2, - • , n are 
transmitted simultaneously from the n transmit antennas. The coefficient a j y is the path 
gain from transmit antenna i to receive antenna j . The path gains are modeled as 
samples of independent complex Gaussian random variables with variance 0.5 per real 
dimension. The wireless channel is assumed to be quasi-static so that the path gains are 

30 constant over a frame and only vary (if at all) from one frame to another. 
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At time t the signal r eJ , received at antenna j , is given by 

n 

where 7] t f represents noise samples. The average energy of the symbols transmitted 
from each antenna (c t i ) is normalized to be 1 In , so that the average power of the 

5 received signal at each receive antenna is 1 and the signal to noise ratio is SNR. 

When block encoding is employed, a block of information bits that arrives in p 
time slots is converted to k symbols that are mapped via a space- time block code onto 
constellation points and transmitted over the n transmitting antennas over the p time slots. 
The receiver receives the transmitted signals and, assuming coherent detection, computes 

10 the decision metric 

P m n 

over all codewords 

and decides in favor of Ihe codeword that minimized the equation (2) is sum. 
15 Encoding in Transmitter 10 

A space- time block code is defined, therefore, by a pxn transmission matrix G 

(p being the number of time slots in a block, and n being the number of transmitting 
antennas). In the context of this disclosure, the entries of the matrix G are linear 
combinations of variables x x 9 x l9 -~x k , and their conjugates, such that the matrix G 
20 satisfies the equality 

<?G = W+W+-+M)I. (3) 

where 7 M isthe nxn identity matrix For example, (7 S4 represents an arrangement that 
employs blocks of eight time slots each, utilizes four transmit antennas, and may have to 
fomi defined by: 
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See, for example, US Patent 6,088,408 issued July 11, 2000; and particular, see equations 
(27 and (28), which demoastrate a G matrix with linear combinations. When the 
transmission employs a signal constellation with 2 b elements, the encoding of each block 
5 begins with the collection of kb input bits, whereupon the encoder selects k constellation 
signals s l9 --s k for the block. Setting x f - s i for / = 1,2,- , k in G, we arrive at a matrix 

C with entries that are linear combinations of s l , s 2 , - * • s k and their conjugates. 

Thus, while G contains indeterminates x } ,x 2 ,"-x k , matrix C contains vspecific 
constellation symbols (or their linear combinations), which are transmitted from n 

10 antennas as follows. If c t i represents the element in the t th row and the f column of C, 
the entries c tJ J - 1 ,2, - * n are transmitted simultaneously from transmit antennas 
1 ,2, • • * , n at each time slot t = 1 ,2, ■ * • , /; . Hence, the /'* column of C represents the 
transmitted symbols from the i th antenna and the t th row of C represents the transmitted 
symbols at time slot t . Note that C is basically defined using G, and the orthogonality of 

15 G columns (i.e., the attribute of equation (3)) allows the simple decoding scheme 

disclosed herein. Since p time slots are used to transmit k symbols, one can define the 

rate of the code to be k ! p . For example, the rate for the illustrated G# 4 is 1/2. 

Decoding in Receiver 20 

Next, we consider the decoding algorithm. It can be shown that, because G is 
20 orthogonal, the decision metric in equation (2) is the sum of k components each 

involving only the variable x f , / = 1 ,2, - • • ,k ; i.e., each independently involving one of the 

k transmitted symbols. It can be further shown that each component can be computed 
using only linear processing. See, V. Tarokh, J. Jafarkhani and A. R. Calderbank, 
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"Space-time block codes from orthogonal designs," IEEE Trans. Inform. Theory, Vol. 
45, No. 5, pp. 1456-1467, July 1999. Therefore, this is a very simple decoding strategy 
that provides diversity. 

Differential Encoding 

Employing, for example, G g4 and only one receive antenna, the received signal of 
that one antenna, /; , , which can be noted by r, , is related to the constellation symbols 
s l ,s 2 ,s i ,s 4 by 



r, = 



(ignoring noise terms), or 



r 2 =-oc i s 2 +a 2 s x -o^s 4 +a 4 s i , 
r 3 =-a l s 3 +a 2 s 4 +a 3 s, -a 4 s 2 , 
r 4 =-cc s s 4 -a^s, +a 3 s 2 +a 4 s,, 
r } =a,,s* +a 4 s* 4) 
r 6 =-a/ 2 +ay?j" -a 3 j* +a 4 s* 3 , 

r i =-«1^3 -Hx/a +0C 3 S * - a / 2 > 

r s =-a/ 4 -a/ 3 +0& +a 4 s[, 



(1 W 5 Ws ) = (¥ 2 ¥^ > 



where 



a, a, 



a, 



«4 a, 



a. 



-a, -# 4 a 3 a 2 -a 



-a, 



-a, 



-a 2 a 3 
-a, a! 



i 



<*3 

* 

-05 

a* 



«4 



-a. 



In other words, a vector comprising a the set of received signals of a transmitted block of 
signals, that are processed somewhat (in this case, the last four received signals being 
conjugated), is related to a vector comprising the constellation symbols that were 
encoded with matrix G M (the vector s, s 2 s 3 s 4 ) multiplied by a matrix (Q, ) that is related 
to the transmission medium's coefficients. It is noted that the equation (6) relationship is 
not the only one that is possible and that the following equations also hold 



(5) 



(6) 



(7) 
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It is also noted that if the set of constellation symbols s l9 s 2 ,s 3 , s 4 is denoted by S, one 
can rewrite equation (8) as: 

(-r 2 r,r 4 -r 3 -r 6 V 5 V-/7*) = v 2(^ A 

(-r 4 r 3 -r 2 r t -r 8 V 7 * -r 6 V 5 *) =v 4 (£)a 
where v, (S) = * 2 s 3 jr 4 ) r , v 2 (5) = {s 2 - ^ * 4 -* 3 ) r , v 3 (S) = (s 3 -s 4 - s { s 2 f and 
v 4 (5) = ~>y 2 -si / . It is noted further that the four vectors, vi, v 2 , V3, and v 4 , are 
orthogonal to each other. The number of such vectors is equal to the number of columns 
in the G matrix, which is equal to the number of transmitting antennas. 

The thrust of this disclosure is that, at each block, the symbols that are applied to 
the G matrix are related to the symbols applied to the G matrix in the previous block and 
to the input data. Basically, it is an encoding system that includes feedback around the G 
matrix. 

Assuming that we want to employ k symbols to transmit M-bk bits per block (i.e., 
frame), we define a set <7>that consists of 2 hk unit-length vectors , P 2 • - , P oM where each 
vector P w isa kxl vector of real numbers, P w = (P wl P w2 — ,/^) r . That is, the set <P 

consists of 2 bk constellation points that lie on a ^-dimensional unit circle, vriuch allows 
any arbitrarily chosen one-to-one mapping ft to be defined that maps the bk bits of a 

block onto <P. That is, an input set of kb bits maps to, or specifies, a vector P M . , where the 
subscript w is between 1 and 2 W , inclusively. 

The encoding is done by first calculating a k -dimensional vector of symbols 
S = 0] % '•■J i ) r , employing a chosen mapping. Then, indeterminates x,x 2 --x k in Care 
replaced by symbols 5,,%- s k , to establish the matrix C, which is used for transmission 
in a manner similar to a regular space-time block code. That is what the first block 
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(10) 



transmits. What remains, now, is to calculate S = (s l s 2 -~ s k ) r for subsequent blocks so 
that non- coherent detection is possible. 

If S u is the vector of constellation k points that is used for the w th block of M bits, 
then C(S U ) is a matrix that defines what to transmit from each antenna during the 
5 transmission of tW h block. C,(S.),i =1,V-,h isflie column of C\SJ andit 
contains p symbols which are transmitted from the antenna sequentially (in p time 
slots). Given an input of bk bits in block u, the vector P w is specified by the /3 mapping 
described above, in accordance with the principles of this invention the next set of 
symbols S = (s l s 2 • • • s k f , for block ttf-1, is given by: 

£♦1 = P« M (S> P* v 2 (S„ )+ P w , Vi (S K ) + ...+ P wk H (S„ ) 

10 * 

where the P wl term is the ith element (dimension) of the vector P w (that is the mapped 
input data). 

Diflferential Decoding 

Recalling that the received signal, r t , is related to the transmitted signals by 

15 'i= s X«A.i (11) 

(ignoring the noise), and defining 

J? = ('! r 2"Ww , "'2j). (12) 
which assumes p = 2k - as is the case in G M - by using equation (1 1) one can write 
R in terms of the vector S which is used to define the transmitted signals, as follows 

R = S T Cl(a l ,a 2 ,•••,«„), (13) 
(see equation (4)) or R = S T (A{a ] ,a 2 ,- i o^)\A^ pc\,--,a)) (14) 
where the "|" designates concatenation of adjacent matrices, and A(ct, ,a, • • a n ) is the 
ftXrtmatrix 
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a, a 



a, 



« 2 -a, -a 4 



a, 



(%3 a 4 
a 4 -a 3 



-q -a 2 



(15) 



a 2 -o^ J 

Recalling that S u and S„ +I are used for the u"' and (w + 1)"" blocks of bk bits, 
respectively. Using G SA , for each block of data we receive 8 signals. To simplify the 
notation, we denote Hie received signals corresponding to the u th block by r" , r" , ■ ■ ■ r s " and 
the received signals corresponding to the (u + l)" block by r l v+ V 2 y+1 ,---r 8 " +1 . Let us 

With reference to equations (9) and (13) one can construct the vectors 



r 4 



f s r e f i )> 

u % u $ w He 



J \< ~V r 2 ') ; 4 ; 3 '6 '5 



tf*)> 
/?*), 



(16) 



where r,", r 2 " , • • • r g " are the signals that are received at the u' h block, and also observe that 
K = v ? (^ ) r Q . It can also be shown that 0,0.' = 2]T" =i | f I K and, this allows one to 
10 determine that if the product i? a+] R% * is to be developed, its value is: 



(17) 



1=/ 



Recalling that the four vectors, Vj, V2, V3, and V4, are orthogonal to each other , and how 
is constructed in equation (10), results in 



(18) 



15 



where q is 1, 2, 3, or 4 in title above example. Therefore, one can construct a vector 3? 
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<B=(^X* *,. + A 2 * «:* R^K*) 



w4 



(19) 



( 4 > 

i2 



2ll« f 



Because the elements of <P(ie. , the P w vectors) have equal length, to compute 
P K , the receiver can compute the closet vector of ^Pto Once this vector is computed, 
the inverse mapping of j3 is applied and the transmitted bits are recovered, 

5 From the resemblance of equation (19) to analogous fonnula for maximum ratio 

combining, it can be shown mathematically that the above detection method provides 4- 
level diversity assuming 4 transmit and one receive antennas. However, the following 
physical argument also demonstrates this fact. For the multiplicative coefficient 
]T | GCj f in equation (19) to be small, which corresponds to a failure to receive die sent 

10 signal, all | a. |, / = 1 ,2,3 ,4 have to be small. In other words, all sub-channels from die 
four transmit antennas to the receive antenna must undergo fading. Conversely, one can 
say that the decoder suffers from the detrimental effect of fading only if all of the four 
sub-channels have small path gains. This is equivalent to saying that a 4- level diversity 
has been achieved. 

15 The same procedure can be used for more than one receive antenna. For each 

receive antenna j 5 we compute <X> ' using the method disclosed above for <£, assuming 

as if only receiver antenna / exists. Then the closest vector of <P to Y M <f> J 

J-U /=1 

(multiplied by a normative constant) is computed. Subsequently, the transmitted bits are 
computed by applying the inverse mapping of /? . 

20 The above example demonstrates a situation where n = k . When the number of 

transmit antennas Is less than the number of symbols, n < k , the same approach works; 
however, some of the details are different hi what follows we consider an example 
where three transmit antennas, with the space-time block code being: 
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-X, 
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—J 
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4 

+ 

J 

x 



(20) 



When there is only one receive antenna, m =1 , the received signals are related to the 
constellation symbols s x , s 2 , s 3 ,s 4 by 

T[ = or, .v t + a ! ,.v 2 
r 2 =-a l s 2 +a 2 s 1 -a 3 s A , 
r 3 = -oc,5 3 + a,^ + a 3 s t , 
r 4 = -a,5 4 -a/j +a 3 ^ 2 , 

r 5 = a, s' + <x,s* 2 + cc.sl , 
r 6 --^s\ +a 2 s' ~a 3 s' 4 , 
r n =-#[.?* + 02^4 + 0,54, 
r t = -0^4 -^3 +a^. 

One may rearrange the above equations to arrive at 

(1 r iWs wl ) = ( w^s^a, 



(21) 



(22) 



where 



a, a. 



0 a\ 



a. 



a, 



0 



-a. 
0 

-a. 



0 a 3 a 2 -or, 
-a, -cc 2 a* 0 



-a, 



«3 

0 

* 

-a. 



0 

a* 
-a,* 



(23) 



a 2 -a, 0 

One may simply check that Q for G 83 can be calculated from Q for G u , 
10 equation (7), by setting a 4 =0. Therefore, again for each specific constellation symbols 
S, vectors v, ( 5), v, (5), v 3 ( 5), v 4 (5) can create a basis for the 4-dimensional space of any 
arbitrary 4-dimensional constellation symbols and the same encoding and decoding 
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20 



schemes are applicable. Hie only difference in the final result is that ]T ( 4 _ ( \a t fis 
replaced by f as follows: 

3? - (LftX^i? K) 



( 3 



(24) 



Therefore, a 3-level diversity is achieved using three transmit antennas and G g3 as 
expected. 

We have assumed rate half space-time block codes that can be applied to any 
complex constellation set. If the constellation set is real, rate one space- time block codes 
are available and the same approach works. For example, in the case of p - k = 4 , the 
following space- time block code exists for n = 4 : 



G = 



x, 



-x. 



~*4 
*1 



-x, 



7 

* ) 



(25) 



It can be shown that 



(r, r 2 r 2 r 4 r, r* r*r* ) = (s x s 2 s y s A )Q , 
where Q is defined by equation (7). Similar differential encoding and decoding are 
possible if we use the following vectors for J = 1 ,2,3 ,4and R u+l : 



(26) 



41 
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r 2 / 3 /- 4 '| A 
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=(-< 


M r 4 '3 '2 
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=(-< 


K H M W * 
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= (-r; 


'3 r 2 r l 4 


«» * 
'3 


A 2 





(27) 



and 



^f+l jr+1 -f"+0* -(«+])* 



), 



^k+I^W '2 '3 '4 M '2 7 3 

This results in a Ml- diversity, full-rate, scheme for differential detection 

In consonance with the above developments, a transmitter in accordance with the 
principles disclosed herein comprises a mapping module 12 that gathers blocks of bk bit 
and maps them to vectors P w . Those vectors are applied to processing element 13 that 



(28) 
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develops the vector S u +\ , with the help of vector S u > which is derived from the output of 
processing element 13, delayed by one block in element 14. The output of processing 
element 13 is applied to block 15, which computes the symbol in accordance with a 
selected space-time mapping, such as Gg4, and applies the results to n antennas 11-1 
5 through ll-n. 

At the receiver, for each receive antenna (e.g., antenna 21- 1) the received signals 
of a block are formed into vector R u , for example (e.g., in element 22- 1), and with a delay 
element (e.g., 23- 1) the signals Ru+i and R u are developed Those signals are applied to a 

processor (e.g., 24- 1) where the products R u+X 7?* *, R U+X R^, . . . R U + X K * are computed, 
10 which form the vector & 1 for processor 24- 1 (and, generally, <3? J terms are computed in 
the processors 24, where j—1, 2, . . .m). The different ( Q i terms are formed and summed 

in element 25 to fonn ^ J =j 3? ; and minimum distance decoding results in a detected P w 

vector. Lastly, the detected P w vector is applied to inverse mapping circuit 26, which 
recovers bk bits for the transmitted block. 
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